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Abstract 


f/a  Fuel-air  mass  ratio 


The  design  ol  a  high  speed  airbreathing 
propulsion  system  (SCRAMJET)  may  depend  on 
more  precise  information  on  chemical  kinet¬ 
ics  of  the  combustion  processes  than  is  now 
available,  except  for  the  very  simplest 
reactions.  An  experimental  method  was  de¬ 
vised  whereby  the  macroscopic  chemical  re¬ 
action  history  can  be  traced  from  wall 
pressure  measurements  in  a  high  speed  re¬ 
acting  gas  in  a  tube.  The  chemical  reac¬ 
tion  is  started  by  means  of  a  normal-shock 
compression  at  the  tube  entrance  through 
which  premixed  fuel  and  air  is  heated  to 
combustion  temperature.  The  experiments 
were  conducted  with  H2  and  vitiated  air, 
and  the  data  were  compared  with  calculated 
pressure  profiles  using  recent  reaction 
kinetic  information  obtained  from  shock 
tubes.  Good  agreement  between  experiment 
and  calculation  was  found  when  the  flow 
was  steady  and  no  oscillations  were  pres¬ 
ent,  indicating  that  the  kinetics  which 
govern  the  chemical  reactions  in  shock 
tube  experiments  can  be  reproduced  with 
the  standing-wave  experiment.  Wall  effects 
are  evaluated,  In  addition,  a  semi- 
empirical  analysis  was  made  of  the  ignition 
delay  and  recombination,  with  the  assump¬ 
tion  that  the  two-body  reactions  are  es¬ 
sentially  in  equilibrium  during  the  slower 
recombination  period.  The  mathematics  of 
this  approach  is  simple  and  is  potentially 
useful  with  machine  computations  involving 
chemistry  and  mixing,  or  chemistry  and 
aerodynamics.  Finally,  a  short  experiment 
with  methane  combustion  in  the  tube  pro¬ 
duced  similar  phenomena,  indicating  that 
the  method  may  be  useful *f or  determining 
rates  of  heat  release  for  gases  where 
little  is  known  of  the  detailed  kinetics. 

Nomenclature 
A  Cross  section  area 


C 


P 


Specific  heat  at  constant  pressure 


H 

h 

K 

k 

M 

M 

m 

mF 

N 


P 

Q 

R 

Re 

r 

T 

t 

u 

X 


Heat  added  from  combustion 
Enthalpy/unit  mass 

Chemical  reaction  equilibrium  coeffi¬ 
cient 

Chemical  reaction  rate  coefficient 
(appropriate  units) 

Mach  number 

Third  body  in  chemical  reactions 
Molecular  weight 
Mol  fraction 

Moles  of  gas  per  gram  mixture  for 
specific  partial  equilibrium  state 

Moles  of  gas  per  gram  mixture  at 
completion 

Moles  of  gas  per  gram  mixture  at  the 
start  of  chemical  reaction  (at  the 
shock) 

Pressure 

Heat  transferred  to  the  stream 

Universal  gas  constant 

Reynolds  number 

Tube  radius 

Temperature 

Time 

Velocity 

Distance  in  the  tube  along  the  direc¬ 
tion  of  flow 


D  Tube  diameter 

F  Pressure  correction  factor  (Eq.  30) 

f  Coefficient  of  friction 


9  Equivalence  ratio  of  vitiated  air 
supply 

9.  Equivalence  ratio  of  unreacted  fuel  at 
the  start  of  the  chemical  reaction 
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£  Extent  of  recombination,  where 

0  <  4  <  1  as  defined  by  Eq.  (9) 

y  Specific  heat  ratio 
p  Density 

t*  Boundary  layer  displacement  thickness 

p,  Viscosity 

Subscripts 

a  Refers  to  air 

corap  For  compressible  flow 

f  Refer  to  fuel 

inc  For  incompressible  flow 

i,j  Refers  to  chemical  species 

mix  Refers  to  fuel  and  air  mixture 

t  At  total  or  stagnation  condition 

w  At  the  wall 

x  At  a  station  x 

1,2  With  reference  to  the  shock  wave  - 

1  In  front  of  the  shock 

2  Downstream  Of  the  shock 

CO  Free  stream 

Introduction 

The  hypersonic  ramjet  with  its  need  for 
supersonic  combustion  pronounces  a  need  for 
extensive  research  to  determine,  much  more 
precisely  than  heretofore,  the  details  of 
the  chemistry  of  combustion,  For  low  alti¬ 
tudes  where  the  combustor  pressure  is  high 
and  the  chemical  reaction  fast,  the  reaction 
distance  may  be  controlled  by  the  mixing 
processes.  However,  as  flight  altitude  in¬ 
creases,  combustor  pressure  recedes,  and 
detailed  knowledge  of  the  chemistry  of  com¬ 
bustion  is  needed,  since  the  chemical  re¬ 
action  length  may  decide  the  engine  length 
(Ref.  1). 

There  are  other  needs  in  high  speed  com¬ 
bustion;  for  example,  a  simplified  mathe¬ 
matical  description  of  the  chemical  reaction 
could  be  used  to  design  a  nozzle  containing 
a  reacting  gas  and  using  a  method  of  char¬ 
acteristics  with  detailed  chemistry,  or 
analyze  a  mixing  system  in  which  the  gas  is 
reacting  chemically  at  the  same  time.  In 
addition,  a  means  of  studying  overall  re¬ 
action  kinetics  of  other  combustion  systems 
than  H2~air  is  needed,  where  detailed  ki¬ 
netics  are  not  currently  understood  and  are 
considerably  more  complicated.  If  catalysts 
are  added  to  extend  the  range  of  usable  com¬ 
bustion  for  propulsion,  the  chemistry  be¬ 
comes  even  more  complicated  and  incompre¬ 
hensible  . 

In  previous  work  at  the  Arnold  Center 
(Ref.  2),  an  experimental  method  was 


demonstrated  whereby  a  standing-wave  nor¬ 
mal  shock  was  used  to  ignite  an  H2-air  mix¬ 
ture,  and  the  resulting  pressure  change  in 
a  constant  area  tube  was  recorded  from 
pressure  taps  in  the  tube  and  indicated  the 
macroscopic  combustion  history  in  the  flow¬ 
ing  reacting  gas.  In  the  present  work, 

(l)  these  experimental  pressure  profiles 
are  compared  with  pressure  profiles  com¬ 
puted  from  recently  proposed  Ho-air  re¬ 
action  schemes  and  rates,  and  (2)  a  simpli¬ 
fied  mathematical  model  of  the  chemical  re¬ 
action  is  developed,  based  on  the  assump¬ 
tion  that  partial  equilibrium  of  the  chemi¬ 
cal  processes  exists  during  the  recombina¬ 
tion  reactions. 

Theory 

Theory  Of  the  Experimental  Model 

The  purpose  of  the  experiment  is  to  pro¬ 
duce  a  combustible  mixture  and  suddenly 
heat  it  to  ignition  temperature  by  means 
of  a  steady  standing-wave  normal  shock. 

Thus,  the  effect  of  chemical  reaction  in 
the  flowing  stream  can  be  observed  from 
steady-state  wall  static  pressure  measure¬ 
ments. 

A  standing-wave  normal  shock  was  first 
used  by  Nicholls  (Ref.  3),  later  by  Rhodes 
et  al .  (Ref .  4) ,  and  recently  by  Richmond 
(Ref .  5)  as  a  means  of  observing  ignition 
delay.  Reasonable  agreement  was  obtained 
with  other  ignition  delay  measurements, 
although  a  tube  was  not  used  to  confine  the 
flow  downstream  of  the  shock.  In  the 
present  experiments,  the  normal  shock  Is 
established  at  the  entrance  of  a  tube,  such 
that  ignition  temperature  is  produced  down¬ 
stream  of  the  shock  and  the  flow  is  con¬ 
fined  to  a  constant  area  (Fig.  la) .  The 
chemical  reaction  may  be  considered  to 
occur  in  one-dimensional  stream  tubes,  ir¬ 
respective  of  mixing,  if  the  fuel  profile 
and  incoming  gas  properties  are  uniform,  and 
wall  friction  and  heat  transfer  can  be  ne¬ 
glected.  The  effect  of  these  factors  on 
the  experiment  will  be  found  in  the  section 
on  discussion  of  results.  ’This  experimental 
model  is  very  similar  to  the  detonation 
wave  in  a  tube,  as  shown  in  Fig.  lb.  The 
principal  phenomenological  difference  is  in 
the  boundary  layer,  in  which  the  velocity 
is  accelerating  near  the  wall  with  respect 
to  the  wave  in  the  shock  tube,  and  decelerat¬ 
ing  with  respect  to  the  wave  in  the  stand¬ 
ing-wave  experiment. 

Whether  the  combustion  zone  is  separated 
from  the  shock  in  a  detonation  wave  has  been 
debated  (Refs.  6  and  7).  However,  the  ex¬ 
periments  at  the  Arnold  Center  with  various 
methods  of  shock-induced  combustion  have 
shown  a  definite  ignition  delay  zone  in 
which  the  pressure  measured  just  downstream 
of  the  shock  corresponds  with  that  which 
would  exist  if  there  were  no  combustion 
(Refs.  2  and  8).  Hence,  the  present  experi¬ 
ment  produces  a  phenomenon  that  agrees  with 
the  Zel  'dovich-von  Neumann-DSring  theory, 
where  the  shock  wave  produces  an  adiabatic 
pressure  rise  followed  by  heat  generated 
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from  combustion.  This  process  is  shown  in 
Fig.  2,  where  the  peak  compression  point  is 
denoted  as  the  von  Neumann  spike  (pt.  B) . 
constant  area  heat  addition  occurs  on  a 
Rayleigh  line  (B  to  C) ,  and  if  sufficient 
heat  is  added  under  steady  conditions,  the 
thermal  choking  or  Chapman-Jouguet  point  is 
reached . 


a . 


Choke 
Point 

Heat 

-x — — — — -  Transfer 

s  Reaction  Zone 

Standing  Wave  Model 


b. 


Velocity  ot  Wall 
Relative  1o  Shock 
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'l  i'.T’  . .  V:  - 

-Boundary  layer 
Shock  Tube  Detonation  Wave  Model 


Fig.  1  Comparison  of  Shock  Tube  Detonation 
Wave  Model  with  Standing  Wave  Model, 
Referenced  to  the  Shock  Wave 
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Fig.  2  Hugoniot  Curve  of  Detonation  Wave 

and  Nornal-Shock-Induced  Combustion 


Maintaining  the  shock  exactly  on  the 
model  lip  is  physically  very  difficult  with 
heat  release  alone  because  slight  varia¬ 
tions  in  fuel  concentration  or  airflow 
properties  will  cause  fluctuations  in  the 
shock,  resulting  in  spillage  or  swallowing. 
Therefore,  the  flow  restricting  cone  plug 
was  used  to  help  maintain  the  shock-on-lip 
position.  Thus  the  choke  point  was  reached 
through  a  combination  of  chemical  heat  ad¬ 
dition  from  point  C  to  point  D  (Fig.  2) 
and  flow  restriction  with  the  plug  to  reach 
the  choking  point  at  the  exit  at  point  E, 
probably  with  some  continued  chemical  re¬ 
action  . 

Theory  of  the  Chemical  Model 


The  behavior  of  an  inviscid  quasi-one- 
dimensional  reactive  flow  may  be  followed 
quite  accurately  by  using  numerical  inte¬ 
gration  techniques  on  the  set  of  equations 
which  express  the  rate  laws  for  the  chemi¬ 
cal  transformations  and  the  appropriate 
conservation  equations  describing  the  flow. 
Both  the  GASL  (Ref.  9)  and  APL  (Ref.  10) 
computer  programs  are  examples  of  such  an 
approach,  and  indeed  comparisons  made  at 
the  Arnold  Center  (Ref.  11)  show  that  with 
identical  kinetics  these  two  programs 
yield  essentially  equivalent  results.  The 
difficulties  with  these  programs  are  two¬ 
fold.  First,  in  order  for  this  approach  to 
yield  accurate  predictions  for  even  the 
simplest  flow  situations  one  must  know 
which  reactions  are  important  and  one  must 
have  reasonably  accurate  rates  for  all  the 
participating  reactions.  Second,  this 
technique,  while  adequate  for  quasi-one- 
dimensional  flows,  is  prohibitively  time 
consuming  in  other  more  complicated  situ¬ 
ations:  i.e.,  in  a  method  of  characteris¬ 
tics  program  or  in  a  program  which  allows 
diffusive  or  turbulent  mixing.  In  general, 
therefore,  one  must  first  explicitly  prove 
the  applicability  of  the  kinetics  used  in 
any  prograai  of  this  type  and  then  if  possi¬ 
ble  find  appropriate  simplifications  to  the 
kinetic  scheme  (see  lor  example  Snow,  Ref. 
12.  on  the  use  of  the  steady-state  approxi¬ 
mation)  to  make  the  program  tractable  for 
the  more  general  type  of  calculation. 

H2  Kinetics 


In  the  case  of  the  H2-air  reaction,  vari¬ 
ous  investigators  have  selected  chemical 
reaction  schemes  and  rates  that  appeared 
most  probable  based  on  the  then  current  ex¬ 
perimental  evidence.  Examples  are  the  APL 
and  GASL  reactions  and  rates.  Reaction 
histories  computed  from  these  two  programs 
agree  very  well  with  each  other  and  with 
the  experimental  data  for  the  ignition  de¬ 
lay  period  (Ref.  8).  However,  large  dif¬ 
ferences  appear  between  the  APL  and  GASL 
predicted  recombination,  (i.e.,  heat  release) 
histories,  and  in  their  comparison  to  ex¬ 
periment  (Ref.  2). 

Recently,  Getzinger  (Ref.  13)  has  re¬ 
ported  new  rate  constants  for  the  recombi¬ 
nation  reactions  and,  in  particular,  has 
stated  that  only  three  reactions  are 
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important  during  recombination.  Therefore, 
the  APL  computer  program  was  modified  to 
include  only  the  following  reactions: 

1.  The  Shuffle  reactions  (Ref.  10): 


and 


k- 5  =  k5/(Kp5  RT) 


k-6  k6/^Kp6  RT'1 


1.  Ho0  +  H  r  H-  +  OH 


kl  -  1.0  x  10 15  exp  (-2SOOO/RT) 
cc/mole  sec 

k2 

2.  H00  •*-  0  r  OH  +  OH 

k-2 

k2  =  5.0  x  1014  exp  (-18000/RT) 
cc/mole  sec 


where  the  units  on  are  atmosphere,  R  has 
units  of  cc  atm/mole  °K,  and  T  is  in  °K . 
Reaction  (7)  presents  a  problem  in  that  it 
produces  the  new  specie  HOn,  which  must 
either  be  included  in  the  full  scheme  or 
handled  in  some  other  manner.  The  con¬ 
sensus  at  present  (Ref.  13)  is  that,  at 
high  temperatures,  this  specie  is  quickly 
destroyed  by  reaction  with  reactive  radi¬ 
cals.  Therefore,  HO2  was  not  included  as  a 
specie  in  the  kinetic  scheme  but  was  assumed 
to  react  rapidly  with  hydrogen  atoms  to 
produce  two  OH  radicals  through  the  process: 


k3 

3.  H2  +  0  “  OH  +  H 

k-3 

k3  -  7.0  x  1012  exp  (-8500/RT) 
cc/mole  sec  and 

k4 

4.  0„  +  H  r  OH  +  0 

-4 

k4  -  5.0  x  1014  exp  ( -18000/RT) 
cc/mole  sec 

where  R  has  units  of  cal/mole  °K,  T  is  in 
°K,  and  the  reverse  rates  of  these  reac¬ 
tions  are  assumed  to  be  related  to  the  for¬ 
ward  rates  through  use  of  the  equilibrium 
constants, 

k  =  k  /KD 
-m  m  Pm 

2.  The  recombination  reactions  (Ref.  13): 
k5 

5.  H  +  H  +  M  2  H„  +  M 

k-5 

k*.  =  4.0  x  lO1,4  [l  +  4(mFH  0  +  mFjj  )] 

2  2  2  2 
cc/mole  sec 

k6 

8.  H  +  OH  +  M  =  Hn0  +  M 
k-6  2 

15  2  2 

kg  =  5.5  x  10  cc  /mole  sec  and 


7.  H  +  0g  +  M  — ■  HO 2  +  M 

k7  =  1.42  x  1015[l-4mF.,  +  29ir.Fu  _,] 

I  U2  HpU 

2  2 
cc  /mole  sec 


In  the  scheme  that  was  employed,  reac¬ 
tions  (5)  and  (6)  are  assumed  to  be  revers¬ 
ible  with  the  reverse  rates  given  by  the 
expressions , 


8. 


H  ■+■  H02  —  20H 


Therefore,  reaction  (7)  is  assumed  to  be 
irreversible  and  followed  rapidly  by  reac¬ 
tion  (8)  to  produce  the  overall  stoichio¬ 
metric  change: 

2H  +  Og  —  20H 

It  should  be  added  that  initiation  re¬ 
actions  such  as  02  +  M  — •  20  +  M  were  not 
included  in  this  scheme  although  they  are 
undoubtedly  important  to  the  early  process¬ 
es  occurring  in  the  induction  zone  (Ref.  14); 
all  the  comparisons  reported  herein  were 
made  for  shock-induced  combustion  in  viti¬ 
ated  air  which  already  contained  radicals 
(i.e.,  OH  levels  near  10_®  mols/mol  mixture). 

Numerical  integration  was  started  just 
downstream  of  the  normal  shock  for  all  the 
calculations.  The  method  of  determining 
initial  conditions  for  a  specific  calcula¬ 
tion  will  be  described  later  in  this  section. 

The  Partial  Equilibrium  Concept 

Kaskan  and  Schott  (Ref.  15)  and  Schott 
(Ref.  16)  have  discussed  the  requirements 
for  partial  equilibrium  (PE)  in  a  reacting 
system.  In  short,  they  define  a  PE  state 
as  one  in  which  the  available  reversible 
reactions  will  not  allow  full  equilibrium. 
They  also  point  out  that  in  certain  cases 
one  may  treat  a  reacting  system  as  having  a 
partial  equilibrium  state  if  a  number  of 
reversible  reactions,  which  are  in  them¬ 
selves  incapable  of  producing  full  equilib¬ 
rium,  are  rapid  while  all  other  reactions 
are  slow.  Thus,  depending  on  the  differ¬ 
ence  in  the  relative  rate  of  these  two 
groups  of  reactions,  the  entire  kinetic 
system  may  be  more  or  less  accurately 
modeled  by  assuming  that  the  fast  reactions 
define  the  relative  concentrations  of  all 
the  species  and  that  the  slow  reactions 
simply  drive  this  partial  equilibrium  (PE) 
system  through  a  series  of  PE  states  toward 
full  equilibrium. 

In  two  significant  papers,  Schott  and 
Bird  (Ref.  17)  and  Getzinger  ar.d  Schott 
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(Ref.  18)  have  shown  that,  at  least  in  the 
range  of  temperatures  from  2500  to  3200°R , 
the  recombination  kinetics  of  the  hydrogen- 
oxygen  system  may  be  modeled  Quite  adequate¬ 
ly  by  using  the  PE  approach.  In  their  de¬ 
velopment  of  the  concept,  they  assumed  that 
the  induction  period  of  the  reaction  occurs 
without  any  simultaneous  recombination  and 
that  at  the  end  of  the  induction  period, 
and  thereafter,  all  shuffle  reactions  (i.e., 
the  binary  exchange  reactions  1-4)  remain 
in  equilibrium  while  the  recombination  re¬ 
actions  simply  change  the  extent  of  recombi¬ 
nation.  Since  the  shuffle  or  exchange  re¬ 
actions  do  not  cause  a  change  m  the  molec¬ 
ular  weight  of  the  system,  they  define  the 
momentary  PE  state  of  the  system  through 
use  of  a  single  variable  which  represents 
the  instantaneous  value  of  the  average  mo¬ 
lecular  weight  of  the  system.  The  same 
approach  will  be  used  here. 

The  PS  Program 

For  this  program,  the  initial  composi¬ 
tion  of  the  vitiated  mixture  is  defined  in 
terms  of  three  variables:  qcv,  the  vitiated 
stoichiometry;  Of,  the  unreacted  hydrogen 
stoichiometry  at  the  normal  shock  wave:  and 
i;,  the  extent  of  recombination  where 

N  -  N 

C  -  -1 -  (o) 

N  -  N , 


and  Ni  =  moles  per  gram  of  mixture  at  the 
shock,  N  -  moles  per  gram  of  mixture  at  any 
specific  partial  equilibrium  state,  and 
Nf  =  moles  per  gram  of  mixture  if  the  re¬ 
combination  reactions  were  to  go  to  full 
completion  (i.e,,  to  the  point  where  the 
mixture  contains  only  H2O  +  O2  +  inerts  for 
a  lean  mixture,  or  H2O  +  H2  +  inerts  for  a 
rich  mixture) .  It  is  further  noted  that 
the  specification  of  cpv,  pj .  g,  and  T  com¬ 
pletely  describes  the  system  because  the 
partial  equilibrium  composition  is  deter¬ 
mined  by  a  set  of  three  equilibrium  equa¬ 
tions,  which  are  pressure  independent  (since 
they  may  be  represented  entirely  by  bimolec- 
ular  exchange  reactions) .  The  program, 
therefore,  contains  a  subroutine  which  cal¬ 
culates  the  partial  equilibrium  composition 
of  the  mixture  for  the  set  of  variables 
(cpv,  <Pf,  T)  .  For  this  calculation,  one 
must  solve  three  non-linear  equations  in 
three  unknowns.  A  three-variable  h'ewton- 
Raphson  method  was  employed  in  this  sub¬ 
routine  to  perform  the  calculation  iterat¬ 
ively  . 

The  remainder  of  the  calculation  was  per¬ 
formed  in  the  following  manner.  The  experi¬ 
ment  yielded  the  following  properties  of 
the  fuel  and  vitiated  air  mixture  before 
reaction:  total  temperature  (Tt) ,  total 

pressure  (Pt) .  static  pressure  ahead  of  the 
normal  shock  (Pi),  and  static  pressure  im¬ 
mediately  downstream  of  the  lip  shock  (P2) . 
In  addition,  cpf  and  -pv  could  be  determined 
from  experimental  measurements.  Only  two 
of  the  three  experimentally  measured  pres¬ 
sures  ai e  needed  to  perform  a  calculation. 


The  pressures  (Pi  and  P2)  were  chosen,  and 
Ptg  was  calculated  for  comparison  with  the 
experimental  value. 

In  the  program,  an  ideal  gas  assumption 
was  used  with  real  gas  enthalpies.  The  en¬ 
thalpies  of  the  species  were  obtained  from 
five-point  curve  fits  to  the  JANAF  tables 
(Ref.  14),  and  equilibrium  constants  were 
fitted  to  the  equation  (log  Kp  =  A/T  +  B) 
where  A  and  B  are  constants.  A  subroutine 
was  used  to  calculate  the  enthalpy  (H) , 
heat  capacity  ratio  (7),  constant  pressure 
heat  capacity  (Cp) .  and  the  molecular  weight 
(m)  of  the  mixture  at  a  temperature  (T)  each 
time  it  was  required. 

The  flow  conditions  upstream  and  down¬ 
stream  of  the  normal  shock  were  calculated 
in  the  following  manner:  The  stagnation 
enthalpy  (H)  was  determined  for  the  un¬ 
reactive  mixture.  At  this  point,  a  series 
of  temperatures  lower  than  T-j.  were  assumed 
successively  until  the  calculated  Mach  num¬ 
ber  became  greater  than  1.5.  Once  two 
supersonic  input  flow  values  were  obtained, 
a  double-loop  iterative  calculation  was 
used  to  adjust  the  incident  temperature  un¬ 
til  the  calculated  value  of  P2  Pi  for  a 
normal  shock  equaled  the  experimental  value. 
The  shock  calculation  is  described  below. 
Stagnation  pressure  (Pt2^  was  then  calcu¬ 
lated,  assuming  unreactive  isentropic  stag¬ 
nation  for  the  flow  downstream  of  the  shock. 

The  shock  calculation  was  performed  by 
using  the  static  enthalpy,  static  tempera¬ 
ture.  and  molecular  weight  of  the  mixture 
at  two  stations  (one  upstream  and  one  down¬ 
stream  of  the  shock)  to  determine  the  other 
property  changes  across  the  shock.  To  this 
purpose,  one  may  solve  the  equations  of 
conservation  of  mass,  momentum,  and  energy 
for  an  ideal  gas  with  variable  molecular 
weight  and  real  enthalpies  to  directly  yield 
the  density  ratio  between  the  two  stations 
in  terms  of  the  above  variables.  The  re¬ 
sulting  equation  is: 


2RT2-  m2  T2ml 


(10) 


The  incident  velocity 
by  combining  the  mass 
equations  to  yield 


p2 


Pi 


2(h2  -  hL) 

(>Vp i)2  " 


may  then  be  obtained 
and  energy  balance 


1/2 


1 


(11) 
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The  remaining  flow  and  state  variables  may 
be  calculated  with  this  information.  For 
the  first  shock  calculation  discussed,  the 
shock  was  assumed  to  be  unreactive,  and 
therefore  m2  -  .  This  calculation  was 

used  to  set  initial  values  downstream  of 
the  shock  wave  for  both  the  APL  calculation 
with  complete  chemistry  and  the  simplified 
partial  equilibrium  program. 

In  the  PE  program,  the  shock  equations 
were  solved  iteratively  for  a  series  of 
assumed  values  of  4.  incremented  from  0  to 
a  value  near  1.  Thus  all  flow  properties 
were  uniquely  determined  for  specific  values 
of  4.  The  only  remaining  problem  concerns 
the  distance  between  these  stations.  For 
4  =  0  (the  end  of  the  induction  zone) .  the 
approach  of  Rhodes  (Ref.  20)  was  used  as  an 
approximation  to  calculate  the  induction 
delay  time.  For  this  purpose,  the  radical 
mole  fractions  were  assumed  to  be  in  equi¬ 
librium  in  the  vitiated  mixture  at  the  tun¬ 
nel  throat  and  were  used  as  initial  radical 
concentrations  downstream  of  the  shock  with 
the  assumption  that  the  composition  was 
frozen  during  expansion  in  the  nozzle. 

These  concentrations  are  quite  low,  and 
they  were,  therefore,  not  included  in  the 
overall  stoichiometry  of  the  system,  A 
constant  multiplying  factor  was  included  in 
the  equation  for  the  final  concentrations 
so  that  the  length  of  the  induction  zone 
could  be  adjusted  if  desired.  It  was  found 
that  setting  this  factor  at  0.5  yielded  the 
best  correlation  of  all  data.  By  using  the 
induction  time  thus  calculated,  the  length 
of  the  induction  zone  was  found  by  multiply¬ 
ing  it  by  the  average  flow  velocity  in  the 
interval . 

During  recombination,  the  distance  be¬ 
tween  stations  (1)  and  (1  +  1)  having  values 
of  £  which  differ  by  A4  was  calculated  by 
determining  dg/dt  at  each  station,  using 
the  equation 

d4  P2  *  -T,  (~k5mFH  ,  „  ] 

dt  (N.-Nf)m2  H  L  2  SmF0H +  k7mF02  J 

(12) 


The  1/2  in  the  k5  factor  arises  because  the 
concentration  of  hydrogen  atom  changes 
twice  as  rapidly  as  4  in  this  equation  and 
the  original  chemical  notation  is  written 
in  terms  of  H-atom  concentration  change. 
These  values  of  d4/dt  are  then  used  to  de¬ 
termine  the  distance,  from  the  equation 


Ax  = A4(u^ 


(13) 


In  the  equation  for  d£/dt,  it  is  assumed 
that  the  end  state  for  the  recombination 
process  is  a  mixture  in  which  N  =  .  This 

contention  was  checked  in  each  calculation 
by  calculating  the  full  equilibrium  concen¬ 
trations  through  application  of  one  of  the 
recombination  equilibrium  constants.  In  all 


cases  reported  herein,  4  >  0.99  at  equi¬ 
librium.  Therefore,  the  neglect  of  reverse 
reactions  in  the  equation  for  d£/dt  was 
satisfactory  for  comparisons. 

This  PE  approach  may  also  be  used  for 
other  flow  conditions.  In  the  case  of 
constant  pressure  flow,  for  example,  the 
enthalpy  per  gram  of  mixture  remains  con¬ 
stant  and  the  iterative  procedure  involves 
finding  T(4)  such  that  h(4)  -  hinitial-  The 
procedure  is  relatively  simple  because 
through  the  use  of  the  PE  concept  the  prob¬ 
lem  has  been  reduced  to  a  single  reaction 
coordinate  (4)  which  is  independent  of 
pressure.  Thus,  this  technique  is  also 
easily  adaptable  to  a  method-of -character¬ 
istics  or  mixing  analysis. 

Equations  Used  in  the  Boundary 
Layer  Corrections. 

The  equations  of  conservation  of  mass, 
momentum,  and  energy  can  be  used  to  calcu¬ 
late  the  equations  for  flow  and  chemical 
change  in  a  duct  (Ref.  21).  By  differenti¬ 
ating  the  equations,  the  "influence  coeffi¬ 
cient"  may  be  calculated  for  the  effects  of 
variables  on  static  pressure  change  in  the 
flowing  reacting  mixture: 

dP  _  -yM2  dQ  _  yM2  dH  yM2  dA 

P  1  -  M2  CpT  1  -  M2  CpT  1  -  M2  A 

_  TM2Ll  +  (Y  -  DM2]  «.  dx  +  -YM2  dm 

2(1  -  M2)  D  1  -  M2  m 

(14) 

The  molecular  weight  term  may  be  neglected 
since  both  the  calculated  and  experimental 
data  have  molecular  weight  changes  included 
For  undeveloped  flow  in  a  tube,  the  bound¬ 
ary  layer  growth  can  be  considered  in  a 
similar  manner  to  boundary  layer  on  a  flat 
plate  where  boundary  layer  displacement 
thickness  includes  heat  transfer  and  fric¬ 
tional  losses  and  can  be  treated  as  an 
effective  area  change.  Equation  (14)  then 
reduces  to: 

<15) 

P  1  -  M2  CpT  1  -  A 


where  the  dH  term  represents  chemical  heat 
release  as  calculated  and  the  dA  term  the 
effective  area  change  from  friction  and 
heat  transfer.  The  tern,  dA/A  can  be  evalu¬ 
ated  from  (see  appendix) 


where  6*  represents  displacement  thickness 
along  the  wall  of  the  tube  and  the  interior 
flow  may  be  treated  as  an  lsentropic  core 
(Ref.  22).  The  displacement  thickness  can 
be  calculated  for  laminar  or  turbulent 
boundary  layer;  for  example  (Ref.  22), 


6 


(17) 
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Since  transition  is  probable  in  the  experi¬ 
ment,  the  end  effect  should  lie  somewhere 
between  the  laminar  and  turbulent  case, 
although  in  either  case,  the  effect  will  be 
shown  to  be  small. 

The  gas  properties  used  to  determine 
Reynolds  number  were  calculated  from  an  H2- 
air  reaction  kinetic  computer  program, 
where  (1)  the  initial  conditions  at  the 
start  of  the  reaction  were  determined  from 
the  experimental  measurements,  (2)  the 
chemistry  was  confined  to  a  constant  area 
process,  and  (3)  the  reactions  and  rates 
produced  a  computed  pressure  history  which 
was  a  close  approximation  to  the  experi¬ 
mental  pressure  history.  The  tube  was  then 
divided  into  small  increments  (Ax)  such  that 
property  changes  were  small  in  each  incre¬ 
ment,  and  the  wall  viscous  and  heat  transfer 
effects  on  pressure  were  calculated  for  that 
increment  with  the  assumption  that  small 
changes  in  the  increment  were  negligible. 
Differential  equations  were  used  with  the 
assumption  that  small  increments  were  a 
satisfactory  substitution  for  the  differ¬ 
ential  term. 

Three  experimental  boundary  layer  cases 
were  considered: 

1.  No  combustion  and  no  heat  transfer. 

2.  No  combustion  but  some  heat  transfer. 

3.  Both  combustion  and  heat  transfer. 

The  pressure  change  measured  for  the 
cases  1  and  2  where  no  fuel  was  injected 
was  compared  with  the  calculated  values. 
Agreement  between  the  experimental  and  cal¬ 
culated  values  was  used  as  justification 
for  applying  similar  methods  to  the  com¬ 
bustion  with  heat  transfer  data  (case  3) 
in  order  to  evaluate  the  extent  of  these 
effects.  As  will  be  discussed,  the  total 
effect  was  small  (on  the  order  of  1  percent)  . 

Error  Analysis 


Effect  of  Change  in  Hydrogen  Fuel  Con¬ 
centration  on  Temperature  of  the  Mixture 
of  Fuel  and  Air 


From  conservation  of  energy  and  the  as¬ 
sumption  that  mean  specific  heat  values  are 
applicable. 


T.  +T  ( f/ a)  (C  /C  ) 

a  _ Pf  pa 

1  +  (f/a)  (C  ,  C  ) 
pf  pa 


(18) 


By  taking  a  first  derivative  and  using  in¬ 
crement  values  to  approximate  the  differen¬ 
tial. 


AT*m 


cPf  <Ttf  Tta) 

cpa  11  + 


A(f/a) 


(19) 


For  typical  experimental  values  of 
f/a  =  0.0016,  A(f'a)  =  0.00016  (10%  error), 
CPf/CPa  s  13>  and  (Ttf  "  Tta>  =  -1000°R, 

which  yields  ATt  -  =  -2°R.  This  is  not 

lmix 

large  enough  to  affect  the  data  or  calcula¬ 
tions  by  a  measurable  amount.  For  a  30- 
perccnt  change  in  f'a  ratio,  AT-j  ^  =  -6°R. 
When  the  fuel  concentration  is  about  0.3 
equivalence  ratio  (E.R.),  a  10-percent  fuel 
concentration  error  will  cause  a  13°R  tem¬ 
perature  error . 

Pressure  Measurement. 

From  the  approach  discussed  in  Ref.  8, 
it  is  found  that  ignition  delay  is  inverse¬ 
ly  proportional  to  pressure.  Hence  a  1- 
percent  error  in  pressure  will  result  in  a 
1-percent  change  in  ignition  delay.  For 
recombination  time  Ref,  9  indicates  it  is 
proportional  to  p~l-7.  Therefore,  a  1- 
percent  error  in  pressure  will  result  in 
approximately  2-percent  change  in  recombi¬ 
nation  time.  Since  the  pressure  measure¬ 
ment  for  most  cases  is  within  1  percent, 
the  effect  of  pressure  error  on  the  chemis¬ 
try  for  the  experimental  cases  discussed 
can  be  considered  negligible. 

Errors  in  Heat  Transfer 

and  Viscous  Effects 

The  wall  effects  for  the  undeveloped 
flow  in  the  tube  is  very  small  so  that  er¬ 
rors  in  calculation  or  measurements  are 
second-order  perturbations.  However,  at 
some  conditions  they  can  be  of  consequence. 
From  Eq.  (15) ,  the  effect  of  error  in 
ary  layer  thickness  is  represented  by 
equation 

AP  ^  yH  AA  2yM  (AS*) 

P  1  -  M2  A  1  -  M2  (r  -  23*) 

where  6*  is  determined  for  laminar  or 
bulent  flow.  For  a  typical  case  near 
front  of  the  tube,  K  =  0.45.  r  -  0.45 
<5*  —  0.005  in.,  and  y  =  1.3,  so  that 


-  3.5^5*) 

x=o 


AP  " 
p  _ 


bound- 

the 


(20) 


tur- 
the 
in  . , 


(21) 


For  an  error  of  10  percent  in  6*  (AS*  =  0 . 0005), 
then  AP/P  s  0,18  percent,  which  can  be  ne¬ 
glected,  However,  as  beat  is  added  and 
Mach  number  increases  to  0.8  near  the  end 
of  the  tube,  with  5*  also  increasing  to 
0.010,  then  AP/P  s  12(A5*) ,  and  a  10-percent 
will  result  in  a  1.2-percent  error  in  AP'P 
or  roughly  of  the  same  magnitude  as  from  6*. 
From  this  analysis,  it  is  concluded  that 
errors  in  determining  the  wall  boundary 
layer  eifects  on  pressure  profile  will  be  of 
negligible  consequence  for  the  fore  parts  of 
the  tube,  including  the  ignition  delay  and 
early  recombination  processes;  however,  as 
heat  is  added  and  the  flow  approaches  M=  1. 
the  wall  effect  errors  may  he  of  the  same 
order  or  greater  than  that  for  the  effect 
Of  the  entire  boundary  layer  displacement 
thickness  . 
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Equipment  and  Procedures 


The  experiment  consisted  of  heating  air 
in  a  supersonic  wind  tunnel  plenum  to  the 
ignition  temperature  of  the  fuel,  or  higher 
(Fig.  3) .  Heating  was  accomplished  by 
means  of  a  heat  exchanger  to  about  1500°R 
and  to  higher  temperature  by  burning  hydro¬ 
gen  ,  Fuel  was  introduced  into  the  viti¬ 
ated  air  in  either  the  supersonic  portion 
of  the  tunnel  or  in  the  subsonic  section 
just  before  the  tunnel  throat  since  in  this 
section  the  static  temperature  was  less 
than  ignition  temperature  (Fig.  4).  The 
possibility  of  preignition  was  greater  when 
fuel  was  injected  at  the  upstream  point, 
but  longer  mixing  lengths  were  obtained. 
Ignition  of  the  fuel  at  the  injector  creates 
distinctive  phenomena  in  the  tunnel  which 
are  easily  recognized.  Emissions  are  ob¬ 
served  upstream  of  the  shock,  and  shock 
position  is  affected  because  tunnel  Mach 
number  is  reduced.  In  addition,  the  com¬ 
bustion  phenomena  in  the  tube  would  not  be 
observed,  except  as  reignition  of  partially 
quenched  reactions.  These  phenomena  were 
not  observed  in  the  present  experiments; 
hence  preignition  was  not  considered. 


Fig,  3  Supersonic  Combustion  Tunnel  with 
Standing-Vave  Normal  Shock 
Combustion  Experiment 
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Fig.  4  Fuel  Injector  Configurations  Used  to 
Premix  Fuel  Entering  the  Combustion 
Tube 


DISTANCE,  IN. 

Fig.  5  Typical  Fuel  Concentration  Distri¬ 
bution  from  a  Multi-Point  Fuel 
Injector  witli  II2  Fuel.  Area 
Weighted  Average  Equivalence  Ratio 
“  0,72;  Center  Point  E.R.  a  0.65. 
Contour  Lines  Are  lor  E.R.  Values. 

Several  fuel  injector  designs  were  used 
in  an  attempt  to  obtain  a  fuel  concentra¬ 
tion  uniformity  of  ±10  percent,  but  this 
effort  was  not  very  successful.  A  fuel 
distribution  profile,  using  an  area  probing 
technique  from  a  single  point  injector,  was 
used  to  design  a  multiple  point  injector 
which  was  calculated  to  produce  a  nearly 
uniform  profile.  One  of  the  better  fuel 
distributions  is  shown  in  Fig.  5.  It  was 
concluded  that,  although  the  fuel  profile 
from  a  single  injection  point  can  be  meas¬ 
ured,  aerodynamic  interference  and  construe 
tion  problems  arise  in  the  process  of  de¬ 
signing  multiple  injection  points  to  obtain 
a  desired  fuel  pattern.  Not  the  least  of 
these  is  the  drilling  of  small  holes  with 
predictable  orifice  coefficients  that  will 
inject  the  gas  in  self  parallel  streams. 
Alignment  procedures  were  used  before  test¬ 
ing  to  check  the  orifice  flow  direction, 
using  water  injection,  and  to  check  the 
flow  direction  after  installing  the  injec¬ 
tors  in  the  tunnel,  using  air  or  nitrogen. 
These  procedures  allowed  fuel  profile  uni¬ 
formity  to  be  improved.  However,  no  design 
manufacture,  or  installation  procedure  ap¬ 
peared  to  provide  the  desired  fuel  profile 
uniformity.  Hence  a  portion  of  the  results 
is  concerned  with  analyzing  the  effects  of 
nonuniform  fuel  profile. 

The  fuel  mixture  subsequently  passed 
through  the  shock  which  was  established  at 
the  combustion  tube  inlet,  and  at  this 
point  the  temperature  was  increased  to  near 
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stagnation  temperature  by  the  shock  com¬ 
pression.  The  shock  therefore  constituted 
the  ignition  plane  and  a  zero  reference 
point  for  the  reaction  kinetics.  The  flow 
can  then  be  considered  one  dimensional 
downstream  of  the  shock  if  there  is  no 
spillover  and  mass  flow  is  conserved.  If 
chemical  reactions  are  proceeding  in  the 
flow,  the  reactions  can  be  considered  as 
events  in  individual  stream  tubes,  inde¬ 
pendent  of  the  other  stream  tubes.  (Trans¬ 
lational  temperature  overshoot  caused  by 
vibrational  relaxation  was  not  considered 
in  the  analysis  since  the  effects  of  this 
phenomenon  on  chemical  reaction  has  not 
been  well  established  and  none  of  the  ex¬ 
perimental  data  show  that  it  was  a  pressure 
measurable  effect.)  The  changes  in  the  gas 
properties  were  determined  from  wall  static 
pressure  measurements,  and  recorded  on 
manometer  tubes. 


Two  combustion  tubes  were  used  (Jig.  6). 
The  thin-walled  spray-cooled  tube  (Fig.  6a) 
was  used  first.  Later  a  jacketed  water- 
cooled  tube  was  built,  which  included  wall 
temperature  thermocouples  and  cooling  water 
temperature  measurement  (Fig.  6b) .  The 
large  base  area  in  this  tube  was  the  prin¬ 
cipal  aerodynamic  difference. 


PRESSURE  TAPS— \  CONE  PLUG 


I 


a.  Spray-Cooled  Combustion  Tube 


COOLANT 

^PRESSURE  PROBE  i  1 

\  rWALL  THERMOCOUPLE 

v  r  - 
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b.  Jacketed  Combustion  Tube 


To  obtain  data  with  minimum  spillover, 
the  data  were  recorded  for  steady  fuel  and 
air  flow,  the  only  variable  being  the  cone 
plug  position.  As  the  cone  plug  was  re¬ 
tracted,  data  points  were  recorded  until  the 
shock  could  not  be  maintained  on  the  lip. 

The  point  recorded  just  prior  to  shock 
swallowing  was  considered  the  best  one. 
Schlieren  photos  were  taken  to  verify  the 
shock  position. 

Measurements 


Measurements  recorded  consisted  of: 

Fuel  Concentration  Measurement 

These  were  made  by  either  (a)  sampling  the 
mixture  in  the  region  ahead  of  the  combustion 
tube  at  several  points  in  a  plane  perpen¬ 
dicular  to  the  ilow,  and  using  a  thermal 
conductivity  gas  analyzer  to  determine  the 
unreacted  hydrogen  content,  or  (b)  measuring 
total  temperature  profile  in  the  heated  air 
stream  where  cold  hydrogen  had  been  injected, 
so  that  hydrogen  concentration  was  calcu¬ 
lated  as  a  function  of  its  cooling  effect. 

For  the  recording  of  combustion  tube  data, 
the  area  survey'  probe  was  removed,  and  a 
vertical  traversing  probe  was  used  which 
measured  total  temperature,  total  pressure, 
and  sampled  the  gas  on  a  vertical  centerline 
in  front  of  the  combustion  tube. 

Tunnel  Flow  Properties  Measurements 

These  consisted  of  (a)  pressure  measure¬ 
ments  to  establish  the  tunnel  operating  con¬ 
ditions  and  (b)  temperature  recording  to 
determine  the  stagnation  temperature  from  the 
preheater  and  the  amount  of  vitiation  in  the 
combustion  preheater.  A  comparison  of 
thermocouple  data  in  the  test  section  and  in 
the  straightening  vanes  showed  good  agree¬ 
ment;  the  straightening  vane  measurements 
were  thereafter  used.  To  calculate  the  con¬ 
centration  of  free  radicals  in  the  test 
section,  it  was  assumed  that  the  vitiated 
air  was  in  equilibrium  at  the  tunnel  throat 
and  that  tile  gas  expanded  with  frozen  com¬ 
position  to  the  test  section  (Refs.  2  and 
8)  . 


Combustion  Tube  Measurements 


Fig.  6  Diagram  of  Combustion  Tubes  <0.91- 
diam  x  7-in.  Long)  Used  ior  Normal- 
Shock  Standing-Vave  Experiments 

The  procedure  for  a  specific  set  of  test 
conditions  consisted  of:  (1)  Adjusting  the 
flow  constricting  exit  cone  so  that  the 
normal  shock  rested  as  close  to  the  lip  as 
possible  without  swallowing,  when  no  fuel 
was  introduced.  This  provided  a  reference 
pressure  profile  which  included  friction 
and  heat  transfer  and  which  would  then  be 
used  to  compare  with  data  containing  fuel 
chemistry,  and  (2)  Injecting  the  desired 
fuel  quantity  and  again  adjusting  the  exit 
cone  so  that  the  normal  shock  rested  as 
close  as  possible  to  the  lip. 


These  consisted  of  23  wall  static  pres¬ 
sure  measurements  in  the  tube  (a)  with 
swallowed  shock  to  determine  stream  static 
pressure  (using  an  average  from  several 
taps) ,  (b)  with  shock-on-lip  and  no  fuel  to 

determine  the  static  pressure  behind  a  nor¬ 
mal  shock,  and  (c)  with  fuel  and  shoex-on- 
lip  to  determine  the  change  in  pressure 
profile  from  chemical  reaction.  In  addition, 
an  actuated  probe  passing  through  the  cone 
plug  was  used  for  some  experiments  to  meas¬ 
ure  total  pressure  gradient  and  sample  gases 
along  the  flow  direction.  Because  of  dif¬ 
ficulties  with  the  probe  it  was  removed. 

For  other  experiments,  the  core  plug  con¬ 
tained  a  sampling  orifice  in  its  tip.  Equiv¬ 
alence  ratio  was  determined  from  gas  sample 
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measurements  from  the  cone  plug.  Where  the 
centerpoint  fuel  concentration  was  not 
representative  of  the  average,  an  adjust¬ 
ment  was  made  In  the  calculations. 

Discussion  of  Results 

The  primary  measurement  in  the  experi¬ 
ment  was  the  wall  static  pressure  drop 
caused  by  combustion  and  heat  addition.  The 
ignition  delay  and,  in  particular,  the 
shape  of  the  pressure  curve  during  the  re¬ 
combination  processes  were  important.  In 
Figs,  7  and  8,  experimental  data  at  two 
pressure  levels  are  compared  with  the  cal¬ 
culated  curves.  The  calculated  curves  were 
obtained  using  recent  H2-air  reactions  and 
rates  (Ref.  18)  and  two  computation  pro¬ 
cedures:  (1)  a  modified  IBM  7094  program 

based  on  the  work  of  Westenberg  and  Favin 
(designated  APL)  (Ref.  10),  and  (2)  a  sim¬ 
plified  "partial  equilibrium"  program  de¬ 
scribed  in  the  section  on  theory,  As  a 
further  comparison,  older  rates  proposed  in 
Refs.  9  and  10  and  tabulated  in  Table  I  for 
the  same  reactions  were  used  in  the  computer 
program  and  also  plotted  in  Fig  7.  The  ex¬ 
perimental  data  appear  to  follow  the  trends 
Of  the  more  recent  reactions  and  rate  data 
proposed  by  Getzinger  and  Schott.  The  ex¬ 
perimental  data  in  Fig.  9  graphically  show 
the  effect  of  increasing  tunnel  temperature 
on  the  chemical  reaction,  particularly  in 
the  ignition  delay  region. 


Fig,  7  Comparison  of  Various  H2-Air  Pres¬ 
sure  History  Calculations  with 
Experimental  Data  and  the  Partial 
Equilibrium  Approximation,  for  an 
Initial  Pressure  of  9.51  psia. 


Table  I 

Comparison  of  Some  Reactions  and  Rates  for  the  H2-A11'  Reaction 


Chemical  Reaction 

"Forward1  Rate  Equation 

From  Ref  9  <GASl) 

From  Rd  10  <APL)  !  Present  Work 

H  +  02-OH  +0 

3  x  IO14  exp  (-8,  81WT! 

5  x  1014  exp  (-18. 000/RT1 

5  x  10 14  exp  f-lg,  000/RTI 

0  +  H2-0H  +  H 

3  x  10“  exp  (-4  030FT) 

7  x  1012  exp  1-8, 300/RTI 

7  x  1012  exp  f-8,5flO/RT! 

2H  +  M  -  H2  +  M 

1  x  1013  x  A 

2  x  io18  r1 

4x  IO14 

41 VV 

1  +  — - - — 

Fi  u 

H  +  OH  +  M  -  HjO  +  M 

l  x  I016  x  B 

3  x  io19  r1 

5  5  X  ID15 

H  +  O  +  M-OH  +  W 

1  x  1013  x  C 

2x  1018  T'1 

20  +  M  -  02  +  M 

3  x  I014 

2x  IG^r1 

HjO  +  H  -  OH  +  H? 

1  33  x  IQ15  exp  (-10.950/T) 

IO15  exp  (-25,  000/RTI 

1013  exp  (-25,  OOO.'RT) 

+  0  -*  20H 

3  12  x  1015  exp  (-12,5!Om 

5  x  IO14  exp  (-18,  OOOfRT) 

5  x  1014  exp  (-18,  000/RT) 

H  +  Oj  +  ff  -  H02  +  M 

H02  +  H  +  M  -  20H  +  M 

— 2H  +  02  +  M  -  20H  +  M 

1  4?  x  1013 

■ 

„  Alternate  Recombination  Factors  (GASU 

Rft  ■  1.9846  ca(/gm°<  ~ — -  n  ,  1  n  ,  „  „ - 

0  „  .  ,  Factor  Rate  1  Rate  2  Rate  3 

T  *  K  I  ,  _  „ 


Factor 

Rate  1 

Rate  2 

Rate  3 

A 

1 

5 

11  25 

B 

1 

10 

31.7 

C 

l 

10 

31  6 

10 


WALL  STATIC  PRESSURE,  PSIA  WALL  STATIC  PRESSURE,  PSIA 


DISTANCE  FROM  LIP,  IN. 


Fig.  8  Typical  Experimental  Pressure  Profile  for  Ii2~Air  Combustion  in 
a  Tube  When  the  Normal  Shock  Is  Stable,  and  P  «  21.58  psia; 
Comparisons  to  Calculated  Curves  with  Complete  Chemistry  and 
with  the  Partial  Equilibrium  Assumption.  Random  Time  Schlieren 
Photos  Show  Shock  Stability. 


Fig.  9  Combustion  Generated  Pressure  Profiles  in  the  Combustion  Tube, 

Using  H2  Fuel,  and  Showing  the  Effect  of  Increasing  Temperature. 


11 


In  addition,  a  short  experiment  was  made 
with  commercial  methane  (95  percent)  to  de¬ 
termine  if  the  combustion  characteristics 
of  methane  were  such  that  similar  recombi¬ 
nation  phenomena  (i.e.,  static  pressure 
change)  could  be  observed  in  the  tube.  The 
pressure  data  clearly  showed  similar  shaped 
pressure  profiles  to  those  from  H2-air  re¬ 
actions.  Ignition  delay  determined  for  the 
experimental  conditions  (T  =  3100  to  3200°R, 
p  =  0.65  atm,  ti.d.  “  10_4  sec)  was  com¬ 
pared  with  shock  tube  results  from  Glass 
et  al .  (Ref.  23)  and  was  found  to  be  faster 
by  a  factor  of  4  to  5.  Preheater  fuel  re¬ 
quired  to  obtain  these  high  temperatures 
(E.R.  =  0.45)  produced  relatively  high  free 
radical  concentrations,  0.3  percent  by  volume 
(Ref.  24).  These  large  amounts  of  free  radi¬ 
cals  and  impurities  could  possibly  affect 
ignition  delay;  however,  no  attempt  was  made 
to  evaluate  their  effect  at  this  time.  The 
experiment  did  emphasize  the  need  for  unvi¬ 
tiated  air  in  studying  reaction  kinetics  of 
combustion  processes  for  which  there  are 
little  data  on  the  detailed  mechanisms. 

As  stated,  the  primary  purpose  of  the 
work  was  to  (l)  correlate  the  measured 
pressure  history  from  the  chemical  reaction 
in  the  tube  with  the  pressure  predicted  by 
a  chemical  kinetic  model,  and  (2)  find  an 
empirical  or  semi-empirical  kinetic  model 
which  could  be  used  to  describe  the  chemical 
reaction  as  a  close  approximation  to  the 
actual  case  and  yet  be  mathematically  sim¬ 
ple.  Both  of  these  objectives  were  accom¬ 
plished  to  some  degree.  The  limitations  of 
the  experiment, /the  limitations  of  the 
theory,  extent  to  which  the  objectives  were 
accomplished,  and  the  problems  which  need 
solution  for  future  work  in  kinetics  experi¬ 
ments  of  this  type  will  be  examined  in  the 
following  sections. 

The  range  of  conditions  for  the  hydrogen 
reactions  studied  in  the  experiments  in¬ 
cluded  initial  temperature  of  1900  to  2300OR, 
initial  pressure  of  0.6  to  1.5  atm,  and 
fuel  equivalence  ratio  of  0.06  to  0.35. 

Experiment  and  Theory  and 
Their  Limitations 

Heat  Release  Compared  with  Theoretical 

Rayleigh  Line  Heat  Addition 

From  the  Zel  'dovich-von  Neumann-Doring 
theory  of  a  shock  wave  followed  by  heat 
addition,  one  may  calculate  the  gas  flow 
properties  just  downstream  of  the  normal 
shock  from  a  knowledge  of  the  upstream 
pressure,  temperature,  velocity  and  compo¬ 
sition.  For  one-dimensional  constant  area 
flow,  the  quantity  of  heat  that  can  be 
added  to  thermally  choke  the  flow  may  be 
calculated.  The  free-streara  Mach  number  for 
the  tunnel  is  approximately  3.1.  It  follows 
that  the  Mach  number  downstream  of  the  shock 
is  approximately  0.45,  and  a  total  tempera¬ 
ture  rise  of  about  13Q0°R  can  be  reached  to 
thermally  choke  a  flow  when  the  initial 
static  temperature  is  2000°R.  This  tempera¬ 
ture  rise  corresponds  to  an  equivalence 
ratio  for  Hg-air  of  about  0.31,  the  upper 


limit  of  fuel  addition  for  a  completed 
chemical  reaction.  Obviously,  higher  equiv¬ 
alence  ratios  can  be  used  if  the  reaction 
is  not  completed  at  the  end  of  the  tube, 
but,  then,  of  course,  less  information  on 
the  reaction  is  obtained.  About  80  to  90- 
percent  reaction  completion  in  the  experi¬ 
ment  is  considered  good  for  comparison  of 
data  with  theoretical  reaction  kinetic 
schemes  and  rates.  The  amount  of  heat  to 
thermally  choke  may  be  increased  by  raising 
the  tunnel  free-stream  Mach  number.  For  a 
Mach  5  free  stream,  choking  equivalence 
ratio  may  be  increased  to  0.55;  and  a  free- 
stream  Mach  number  of  10  is  needed  to  ac¬ 
commodate  an  E.R.  of  1.0.  In  addition,  a 
larger  amount  of  heat  may  be  added  if  the 
initial  gas  temperature  is  higher,  the  tem¬ 
perature  ratio  to  choke  being  determined 
from  Rayleigh  line  considerations.  An  in¬ 
finite  amount  of  heat  may  be  added  to  a 
constant  pressure  combustion  process  with¬ 
out  causing  aerodynamic  choking  of  the  flow, 
even  though  the  Mach  number  may  decrease  to 
1.0  or  less.  However,  the  flow  area  must 
change  very  precisely  in  order  to  maintain 
the  constant  pressure,  creating  a  difficult 
hardware  problem.  Other  duct  shapes  may  be 
usable  for  studying  reaction  kinetics,  such 
as  the  expanding  conical  area.  The  selec¬ 
tion  of  a  configuration  should  be  made  such 
that  the  reaction  does  not  quench  nor  re¬ 
strict  the  entering  mass  flow. 

Because  of  the  Mach  number  and  thermal 
choking  limitation  of  the  experimental 
model,  the  present  approach  is  to  correlate 
experimental  and  theoretical  data  for  lean 
mixtures  only.  Later  experiments  may  be 
used  to  verify  the  reactions  for  higher 
E.R.'s  with  a  configuration  which  is  not 
equivalence-ratio  limited. 

Wall  Heat  Transfer  and  Friction 

The  effects  of  wall  heat  transfer  and 
friction  were  calculated  by  the  method  dis¬ 
cussed  in  the  theory  section.  Pressure 
data  were  obtained  for  the  case  where  (1) 
no  fuel  was  injected,  (2)  the  shock  was 
maintained  on  the  lip,  and  (3)  inlet  air 
temperature  was  varied.  Pressure  profiles 
were  then  calculated  for  these  cases  for 
both  laminar  and  turbulent  boundary  layers. 

A  plot  of  the  data  and  calculations  is 
shown  in  Fig.  10.  Since  transition  to  tur¬ 
bulent  flow  is  probable,  it  would  be  ex¬ 
pected  that  the  pressures  measured  at  the 
end  of  the  tube  would  lie  between  those 
calculated  for  laminar  and  turbulent  flow. 

The  accuracy  of  the  data  and  the  wall  ef¬ 
fects  calculations  is  not  such  that  laminar 
or  turbulent  boundary  layer  can  be  specified. 
The  calculated  effect  is  approximate  and, 
for  three  typical  experimental  cases,  was 
found  to  be  approximately  1.0  percent  at 
the  end  of  the  tube  for  laminar  boundary 
layer  and  1.5  percent  for  turbulent. 

The  region  of  ignition  deLay  and  initial 
recombination  reactions  is  found  in  the 
first  half  of  the  tube.  The  boundary  layer 
here  is  thin,  and  could  be  considered  negli¬ 
gible  or  could  be  used  to  adjust  the  data. 
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Fig.  10  Measured  and  Calculated  Wall  Viscous  and  Heat  Transfer  Pressure 
Effect  in  the  Combustion  Tube  When  No  Fuel  Was  Present 


In  Fig.  8  is  shown  an  experimental  pressure 
profile  that  has  been  corrected  for  laminar 
boundary  layer  displacement  thickness  and 
wall  heat  transfer.  The  corrected  values 
tend  to  more  closely  approximate  the  slope 
of  the  calculated  pressure  profile. 

Fuel  Concentration  Nonuniformity 

A  typical  fuel  concentration  map  is  shown 
in  Fig.  5.  For  some  data  points,  vertical 
traveling  probe  fuel  profiles  from  traverses 
in  front  of  the  combustion  tube  showed  even 
greater  variation.  It  can  be  shown  (Ref. 

25)  that  the  effect  on  wall  static  pressure 
of  a  one-dimensional  flow  with  poor  heat 
distribution  in  a  tube  of  flowing  gas  is 
negligible,  providing  the  pressure  at  any 
one  station  can  be  assumed  to  be  constant 
across  the  duct.  In  addition,  the  computer 
programs  with  reaction  kinetics  can  be  used 
to  establish  that  the  kinetic  history  is  not 
sensitive  to  fuel  concentration  for  lean 
mixtures.  It  is  suspected  in  the  present 
case  that  constant  pressure  across  the  duct 
was  not  the  case  when  fuel  distribution  was 
extremely  asymmetric,  and  that  this  might 
explain  the  apparent  "slow"  recombination 
reactions  obtained  for  a  few  data  points, 
examples  of  which  are  shown  in  Figs.  11  and 


12a.  A  simplified  calculation  of  an  ex¬ 
treme  case  where  the  fuel  was  assumed  to  be 
concentrated  in  10  percent  of  the  flow  next 
to  the  wall  and  the  pressure  taps  were  lo¬ 
cated  on  the  opposite  wall  indicated  that  a 
static  pressure  maximum  lag  in  the  tube 
might  be  as  great  as  0.2  inch  near  the  start 
of  the  reactions,  with  a  maximum  of  0.4  inch 
as  the  flow  approached  M  -  1.  Since  no  fuel 
profile  measurements  showed  a  maldistribu¬ 
tion  of  this  extremity,  the  maximum  effect 
calculated  is  not  possible  for  steady-state 
conditions.  However,  it  is  possible  that 
nonsteady  fuel  fluctuations  may  have  an  im¬ 
portant  effect  on  the  pressure  profile. 

Expelled  Shock  at  the  Lip  with  Spillover 

As  the  flow  constricting  cone  plug  is 
retracted  during  the  experiment,  the  shock 
moves  closer  to  the  lip.  The  condition  most 
closely  approximating  the  theory  occurs  when 
the  shock  is  very  close  to  the  tube  inlet 
and  spillovex-  is  negligible.  When  appreci¬ 
able  spillover  is  present,  the  data  can  not 
be  used  to  correlate  with  the  theoretical 
one-dimensional  flow  model  with  chemical 
reaction.  Data  of  this  type  are  shown  in 
Fig.  11  and  are  discussed  in  Ref.  2. 


Fig.  11  Typical  Combustion  Tube  Pressure  Profile  with  Expelled 

Normal  Shock,  Compared  to  Calculated  Pressure.  Schliercn 
Photos  Recorded  at  Random  Time. 
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Shock  Oscillation  during  Combustion 


It  was  observed  during  the  experiments 
that  at  some  fuel  flows  the  normal  shock 
did  not  remain  on  the  lip  of  the  tube  but 
oscillated  in  a  fore  and  aft  direction  and 
at  an  unknown  (not  measured)  frequency 
(Figs.  12a,  b,  c) .  The  shock  movement 
varied  from  fluctuations  that  were  hardly 
perceptible  to  large  movements  covering 
nearly  the  full  length  of  the  tube  out  to 
about  1/2  inch  in  front  of  the  lip.  The 
magnitude  of  the  movement  appeared  to  be 
roughly  proportional  to  fuel  flow  for  a 
given  set  of  flow  conditions,  The  oscilla¬ 
tions  were  first  observed  when  the  shock 
suddenly  seemed  to  disappear.  If  the  ma¬ 
nometer  board  showed  that  the  shock  had  not 
swallowed,  it  was  noted  that  the  oscilla¬ 
tion  produced  a  peculiar  pressure  profile 
which  indicated  a  "slow1’  recombination  his¬ 
tory  (Fig.  12a)  or  a  pressure  rise  instead 
of  a  drop  near  the  lip  (Figs.  12b  and  12c). 
Furthermore,  from  successive  spark  schlieren 
photos,  the  shock  could  be  observed  at  vari¬ 
ous  positions  during  the  oscillation  (Figs. 
12a,  b,  c)  . 

The  wall  pressure  data  were  considered 
worthless  for  chemical  kinetic  studies  when 
the  oscillations  were  present.  The  pressure 
profile  shown  in  Figs.  12b  and  c  for  an 
oscillating  shock  condition  can  be  analyzed 
qualitatively  if  it  Is  assumed  that  the 
pressure  measured  on  the  manometer  tubes  is 
a  time-mean-average  of  the  instantaneous 
pressures.  This  is  shown  in  Fig.  13.  When 
the  shock  is  expelled,  the  lip  pressure  is 
high  because  of  subsonic  flow  divergence 
behind  the  normal  shock,  and  less  pressure 
change  from  heat  addition  occurs;  when  the 
shock  is  swallowed,  the  low  free-stream 
static  pressure  is  sensed.  The  manometer 
tube  registers  a  pressure  which  is  some 
function  of  the  duration  of  the  pressure 
exposure  time.  In  this  analysis,  it  was 
assumed  that  the  velocity  of  oscillation 
was  relatively  low  compared  with  the  chemi¬ 
cal  reaction  speed,  so  that  the  effect  of 
the  oscillation  on  chemical  reaction  at  any 
one  stage  of  the  cycle  is  small. 

Several  possible  causes  of  the  oscilla¬ 
tions  were  investigated: 

1.  In  the  first  experiments,  a  thin-wallcd 
spray-cooled  tube  was  used,  as  compared 
with  the  thick-walled  water-cooled  tube 
in  the  later  experiments  (Fig.  6) .  Os¬ 
cillation  did  not  appear  in  the  thm- 
walled  tube  until  the  combustion  heat 
release  approached  that  required  for 
thermal  choking.  An  analysis  of  the 
blunt  base  flow  region  was  made  to  de¬ 
termine  if  the  flow  in  this  region  was 
conducive  to  pressure  fluctuations  which 
might  affect  the  choking  ilow  area,  then 
the  mass  flow,  and  hence  the  amount  of 
flow  spillover  at  the  lube  entrance.* 

This  effect,  appears  to  be  a  possibility, 
but  no  definite  conclusions  could  be 


♦Analysis  made  by  R.  C.  Bauer  and 
J .  Garner . 
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Fig.  13  Qualitative  Analysis  Of  Combustion 
Tube  Pressure  Profile  during  Oscil¬ 
lations 


drawn  from  the  base  flow  analysis.  How¬ 
ever.  since  the  base  area  was  the  princi¬ 
pal  difference  between  the  two  tube  mod¬ 
els,  there  is  a  strong  possibility  that 
the  large  base  area  plus  a  combination 
of  exterior  and  interior  aerodynamics 
were  contributing  factors. 

2.  Air  stream  profile  nonuniformity  was  con¬ 
sidered  but  could  not  be  substantiated 

by  pressure  or  temperature  proiile  data 
m  the  air  stream.  Static  pressure, 
total  pressure,  and  total  temperature 
measurements  indicated  very  little  non- 
uniformity  and  less  than  that  observed 
in  the  thin-walled  tube  experiments. 

3.  Fuel  pro! lie  nonuniformity  was  considered. 
From  the  one-dimensional  flow  analysis 

in  Ref.  25,  it  was  shown  that  nonuniform 
heat  addition  would  produce  the  same 
static  pressure  change  along  the  duct 
axis  it  the  static  pressure  at  any  sta¬ 
tion  were  considered  constant  across 
the  duct.  However,  as  mentioned  previ¬ 
ously,  nonuniform  fuel  distribution  can 
result  m  some  variation  in  pressure 
profile.  It  is  conceivable  that  a  fluc¬ 
tuating  fuel  concentration  might  trigger 
oscillating  pressures.  However,  from 
total  temperature  measurements  in  the 
tunnel  stream,  there  was  no  significant 
change  in  small  temperature  fluctuations 
when  fuel  was  injected  as  compared  with 
the  no-luel  case.  The  thermocouple  and 
potentiometer  sensitivity  was  about  10 
cps,  so  that  only  lower  frequency  fluc¬ 
tuations  could  be  effectively  sensed. 

d .  Boundary  layer  separation  with  resulting 
pressure  rise,  cooling  of  the  gas  near 
the  wall,  and  inhibition  of  the  chemical 
reaction  was  considered  but  could  not  be 
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justified  on  the  grounds  that  the  pres¬ 
sure  gradient  in  the  tube  is  a  favorable 
one  and  not  conducive  to  separation. 

5.  Changes  in  reaction  kinetic  rates  were 
considered  as  a  possible  result  of 
changes  in  the  reaction  environment  but 
could  not  be  justified.  Below  1800OR, 

a  shift  in  the  important  chemical  re¬ 
actions  may  be  expected,  as  will  be  the 
Case  for  very  rich  mixtures  (above  E.R, 

"  4) .  Neither  of  these  conditions  ex¬ 
isted  in  the  experiment. 

6.  It  is  also  possible  that  the  oscilla¬ 
tion  may  be  the  result  of  disturbances 
caused  by  the  heat  generating  chemical 
reactions  which  may  trigger  either 
transverse  or  longitudinal  wave  oscilla¬ 
tions.  These  phenomena  are  discussed 

by  Strehlow  and  Fernandes  (Ref.  26)  for 
transverse  waves  and  Fickett  and  Wood 
(Ref.  27)  for  longitudinal  waves.  Fur¬ 
ther  investigations  would  be  needed  to 
ascertain  the  cause  and  characteristics 
of  the  oscillations. 

Therefore,  the  only  experimental  effect 
to  which  the  ''slow"  recombination  history 
could  be  attributed  was  shock  oscillations 
or  inlet  flow  spillage.  For  one  group  of 
experimental  data  at  the  higher  pressure 
level  (90  psia  in  the  plenum) ,  some  data 
points  showed  good  agreement  with  the  cal¬ 
culated  pressure  profile  and  others  showed 
"slow"  recombination  and  poor  agreement. 

A  careful  examination  of  the  data  and  the 
multiple  photographs  of  the  shock  showed 
that  the  better  correlations  could  be  asso¬ 
ciated  with  little  or  no  shock  movement. 

Limitations  of  the  PE  and  Full 
Numerical  Integration  Models 
for  the  Reaction  Zone 

The  APL  program  with  the  inclusion  of 
recently  determined  kinetic  rates  appears 
to  predict  the  observed  pressure  profiles 
quite  well  over  the  range  of  initial  con¬ 
ditions  studied  in  the  current  experiments. 
In  addition,  it  is  quite  apparent  from  the 
results  that  the  occurrence  of  HO2  in  the 
reaction  scheme  may  be  handled  adequately 
by  the  technique  used  in  this  study  over 
the  range  of  experimental  conditions. 

One  change  was  made  in  the  APL  program. 
As  originally  written,  the  program  forced 
mass  conservation  at  each  step  by  calculat¬ 
ing  the  rate  of  change  of  the  hydrogen  and 
oxygen  atom  concentrations  from  the  net 
rate  Of  change  of  the  other  species.  It 
was  noted  that  at  the  start  of  some  calcu¬ 
lations  there  were  sporadic  changes  in  the 
radical  concentrations.  Since  the  rate  of 
formation  of  radicals  is  proportional  to 
the  radical  concentrations,  this  method  of 
handling  the  residuals  could  have  led  to 
these  calculated  instabilities.  Therefore, 
the  program  was  altered  to  carry  the  resid¬ 
uals  in  the  species  O2  and  H2O .  This 
change  completely  eliminated  the  instabili¬ 
ty. 


The  partial  equilibrium  program  appeared 
to  work  satisfactorily  at  high  temperatures 
and  low  pressures  but  did  not  yield  as  good 
an  agreement  at  low  temperatures  and  high 
pressures.  In  these  lean  mixtures,  the 
rates  of  recombination  reactions  are  pri¬ 
marily  controlled  by  the  hydrogen  atom  con¬ 
centration,  and  this  concentration  was 
found  to  be  very  sensitive  to  the  extent  of 
recombination.  An  examination  of  the  APL 
program  results  showed  that  the  hydrogen 
atom  concentration  goes  through  a  distinct 
maximum  near  the  end  of  the  induction  zone. 
However,  at  high  pressure  and  low  tempera¬ 
ture,  this  maximum  is  considerably  lower 
than  that  predicted  for  £  =  0  using  the 
partial  equilibrium  program.  Values  of 
COmn  are  plotted  in  Fig.  14  as  a  function 
of  the  shock  pressure  and  temperature  (?2, 
T2)  .  The  values  of  (£)j,,ln  for  this  plot 
were  determined  by  first  running  an  APL  and 
PE  program  for  identical  conditions.  The 
output  was  when  examined  to  find  the  maxi¬ 
mum  hydrogen  atom  concentration  predicted 
by  the  APL  program.  This  value  was  then 
used  to  determine  an  equivalent  value  of 
(4)min  in  the  PE  program.  These  are  the 
values  identified  in  Fig.  14  as  ( ^) min •  *n 
addition  to  these  points,  a  reasonable  set 
of  constant  (4)min  lines  are  included  to 
indicate  trends,  As  (£)rain  increases,  the 
agreement  between  the  full  numerical  inte¬ 
gration  scheme  and  the  PE  scheme  becomes 
progressively  poorer,  especially  near  the 
end  of  the  induction  zone.  Typical  plots 
of  H-atom  concentration  versus  time  from 
the  two  programs  are  shown  in  Fig.  15. 


INITIAL  STATIC  TEMPERATURE,  °R 


Fig.  14  Evaluation  of  H-Atoin  Concentration 
as  a  Measure  of  Ignition  Delay, 

£min  Plotted  as  a  Function  of  Ini¬ 
tial  Temperature  and  Pressure  at 
the  Peak  H-Concentration .  High 
Values  of  £nin  Indicate  H-Concen- 
tration  Is  Not  a  Good  Indication  of 
Ignition  Delay. 
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H-ATOM  CONCENTRATION,  MOLS/GM  MIX  H-ATOM  CONCENTRATION,  MOLS/GM  MIX 


Fit?.  15  Comparison  of  H-Atom  Concentration 
for  Complete  Chemical  Calculation 
and  Partial  Equilibrium  Calculation 


The  PE  program  approximation  to  estimate 
the  induction  zone  length  appears  to  be 
reasonably  good  at  high  temperature  and  low 
pressure  but  probably  should  be  modified 
at  low  temperature  and  high  pressure  where 
the  recombination  reactions  occur  to  a  sig¬ 
nificant  extent  before  the  exchange  re¬ 
actions  reach  equilibrium.  This  overlap 
of  the  induction  zone  and  recombination 
zone  represents  a  real  limit  to  the  utility 
of  the  PE  model  at  low  temperature  and  high 
pressure  and  warrants  further  investigation 
if  a  semi -empirical  approach  is  desired  for 
these  conditions, 

Extent  of  Accomplishment 

Recently  determined  reactions  for  the 
H2~air  combustion  and  their  rates  were  used 
to  produce  calculated  theoretical  curves  of 
pressure  history  in  a  constant  area  duct. 

In  general,  the  calculated  pressure  his¬ 
tories  show  better  agreement  with  the  ex¬ 
perimental  data  than  calculations  using 
older  proposed  reactions  and  rates. 

Since  the  new  rates  were  taken  from 
shock  tube  experiments,  the  experiments  and 
analysis  constitute  a  verification  of  the 
data  obtained  from  the  shock  tube  experi¬ 
ments  . 

The  theory  that  partial  equilibrium 
exists  during  recombination  for  reactions 
at  low  pressure  and  the  higher  temperatures 
appears  to  be  a  valid  one.  It  was  applied 
as  a  simplified  kinetic  scheme  consisting 
of  an  ignition  delay  portion  based  on  cor¬ 
relation  of  ignition  delay  data  from  Hg- 
air  computer  programs  (Ref.  20),  and  a 
recombination  portion  based  on  three  re¬ 
combination  reactions.  The  pressure  pro¬ 
file  calculated  from  the  PE  relations 
showed  relatively  close  agreement  with  the 
full-scale  reaction  scheme  computations 
using  complicated  mathematics  and  the  ex¬ 
perimental  data  (Fig.  7) .  The  experiments 
were  limited  to  lean  mixtures  up  to  E.R. 

=  0.35:  however,  theory  and  the  work  of 
others  (Ref.  13)  indicate  that  the  PE  as¬ 
sumption  is  valid  from  E.R.  =  0.17  to  4.0, 
The  graph  in  Fig.  14  indicates  that,  for 
lean  mixtures,  there  is  a  high  pressure  and 
low  temperature  range  that  is  not  applicable 
to  the  PE  approach.  A  different  approach 
may  be  necessary  in  arriving  at  semi- 
empirical  equations  describing  the  recombi¬ 
nation  history  for  other  fuels. 

The  initial  gas  composition  included 
free  radicals  and  products  from  the  vitia¬ 
tion  heating.  Hence  the  current  success 
in  predicting  the  reaction  history  from  a 
knowledge  of  the  initial  composition  and 
free  radicals  provides  a  further  informa¬ 
tion  background  for  predicting  the  effect 
of  contaminants  entering  a  combustion  cham¬ 
ber  while  the  air  is  heated  by  vitiation. 
This  is  of  particular  interest  to  those  in¬ 
terested  in  supersonic  combustion  ramjet 
ground  test  facilities  where  the  means  of 
achieving  flight  air  temperature  in  the 
combustor  may  he  by  a  vitiation  process. 
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Problems  lor  future  Work 

With  the  assumption  that  the  work  accom¬ 
plished  so  far  justifies  improving  the 
method  and  extending  its  application  to 
other  combustion  processes,  these  aspects 
deserve  consideration: 

1.  To  eliminate  the  effect  of  contaminants 
and  free  radicals  on  the  reaction,  at 
least  in  the  initial  stages  of  investi¬ 
gation,  a  nonvitiated  air  supply  should 
be  used.  This  can  be  accomplished  for 
a  fuel  with  an  ignition  temperature 
equal  to  or  less  than  hydrogen  by  using 
an  electric  resistance  heater  (to  2500°R) 
or  a  storage  heater. 

2.  Methods  of  improving  the  uniformity  of 
the  fuel  profile  should  be  developed  so 
that  the  flow  entering  the  shock  is  more 
nearly  one  dimensional. 

3.  Higher  equivalence  ratios  can  be  run 
without  thermal  choking  if  the  Mach  num¬ 
ber  of  the  tunnel  is  increased.  This 
will  cause  the  Mach  number  downstream  of 
the  shock  to  be  lower  so  that  a  larger 
heat  release  will  be  required  to  ther¬ 
mally  choke  the  flow.  Another  method 
would  be  to  let  the  combustion  zone  area 
expand;  however,  the  constant  area 
boundary  condition  is  then  lost.  To  es¬ 
tablish  a  contour  for  another  thermal 
heat  release  process,  i.e.,  constant 
pressure,  would  not  be  impossible,  but 
less  likely  to  be  as  accurate  as  the 
constant  area  process.  The  amount  of 
heat  that  could  be  added  to  a  constant 
pressure  process,  however,  would  not  be 
limited.  Other  methods  would  include 
expanding  the  flow  area  to  prevent  ther¬ 
mal  choking,  or  increased  initial  tem¬ 
perature  where  the  choking  temperature 
ratio  is  fixed  from  Rayleigh  line  cri¬ 
teria  . 

4.  Other  fuels  can  be  investigated  experi¬ 
mentally,  although,  since  much  less  is 
known  about  other  combustion  processes 
than  H2-air,  the  equations  describing 
the  reactions  would  likely  contain  more 
empiricism.  A  short  test  was  made  with 
commercial  (95  percent)  methane-air  com¬ 
bustion  in  the  existing  tunnel,  and  al¬ 
though  a  high  degree  of  vitiation  was 
necessary  to  obtain  combustion  (approxi¬ 
mately  3300°R  tunnel  plenum  temperature), 
the  shapes  of  the  pressure  profiles  were 
similar  to  those  obtained  from  hydrogen 
combustion . 

5.  If  catalysts  are  mixed  with  the  fuel  or 
air,  then  it  would  be  expected  that 
their  effect  on  the  combustion  process 
can  be  determined  by  changes  in  the  pres¬ 
sure  profile  in  the  tube.  Information 

of  this  type  is  of  interest  in  super¬ 
sonic  ramjet  combustion  research,  where 
the  upper  limits  of  flight  altitude  and 
Mach  number  may  be  restricted  because  of 
the  time  required  for  the  combustion. 
Catalysts  which  speed  the  reaction  may 
extend  the  range  of  operation  in  the 
flight  envelope  of  these  engines . 


6.  Because  the  presence  of  flow  oscilla¬ 
tions  is  incompatible  with  the  one¬ 
dimensional  flow  analysis  approach 
toward  reaction  kinetics,  this  effect 
should  be  removed.  Fuel  profile  uni¬ 
formity  may  be  improved  by  using  mix¬ 
ing  theory,  and  the  model  may  be  re¬ 
designed  to  reduce  the  occurrence  of 
oscillations.  However,  transverse  and 
longitudinal  oscillations  observed  in 
traveling  detonation  waves  may  also 
occur  under  similar  circumstances  in 
the  standing  wave,  and  it  may  be  that 
the  oscillations  cannot  be  completely 
eliminated.  Further  studies  of  the 
phenomena  may  reveal  ways  of  controlling 
it . 

Comparison  of  Results  with 
Shock  Tube  Experiments 

From  a  one-dimensional  flow,  ideal  gas, 
point  of  view,  the  detonation  wave  in  a 
tube  and  the  standing-wave  normal-shock 
combustion  are  identical.  From  a  physical 
point  of  view,  there  are  several  differ¬ 
ences.  Since  data  obtained  from  detonation 
experiments  (Ref.  18)  were  used  to  corre¬ 
late  theoretical  calculations  with  the 
standing  wave  experimental  data,  a  compari¬ 
son  of  the  two  phenomena  is  of  interest. 

Heat  Transfer  and  Viscous  Effects 

Because  the  boundary  layer  for  the  deto¬ 
nation  tube  is  accelerating  with  respect  to 
the  wave  (Fig.  1),  the  friction  term  will 
be  positive.  Heat  transfer  will  still  be 
toward  the  wall  if  the  wall  is  cold  (Ref. 
28).  From  Eq .  (14),  it  can  be  seen  that 
the  two  effects  will  be  cumulative,  rather 
than  opposite  and  compensating  as  m  the 
case  of  the  standing  wave,  and  the  dis¬ 
placement  boundary  layer  will  be  negative 
(expanding).  These  wall  effects  will  also 
be  influenced  by  tube  diameter. 

Shock  Wave  Stability 

The  standing-wave  normal-shock  stability 
depends  on  the  steadiness  of  the  entering 
flow  and  on  the  tendency  of  the  shock  wave 
to  oscillate  when  combustion  occurs.  On 
the  other  hand,  the  flow  entering  the  trav¬ 
eling  detonation  wave  does  not  contain  the 
turbulence  that  is  characteristic  of  the 
flowing  gas  system,  and  gas  mixture  uni¬ 
formity  is  easily  obtained.  However, 
transverse  and  longitudinal  oscillations 
are  observed  in  detonation  experiments, 
and  the  one-dimensional  analysis  of  re¬ 
action  kinetics  can  not  then  be  applied 
(Ref,  29).  It  is  conceivable  that  a  tunnel 
and  tube  model  redesign  for  the  standing- 
wave  experiment  could  produce  a  relatively 
uniform  flow  and  fuel  mixture  and  reduce 
the  base  flow  problems,  thus  reducing  the 
probability  of  the  oscillations.  However, 
the  standing-wave  experiment  may  be  sub¬ 
ject  to  similar  oscillation  phenomena 
which  have  not  been  completely  explored  in 
the  present  experiments. 
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Experimental  Measurements 

Microsecond  response  time  is  needed  for 
shock  tube  measurements,  and  ior  this 
reason,  optical  measurements  are  the  most 
reliable.  The  shock  tube  equipment  is 
relatively  simple,  but  the  measurement 
equipment  and  techniques  become  quite  elab¬ 
orate.  On  the  other  hand,  the  standing- 
wave  combustion  tube  measurements  may  be 
made  by  relatively  simple  instruments,  such 
as  the  manometer  board;  but  the  steady  flow 
requirements  create  the  need  for  relatively 
complicated  equipment,  such  as  large  air 
and  fuel  supplies,  air  heaters,  and  controls. 

Overall  Kinetic  Consid erations 

At  the  present  stage  of  development  of 
the  standing-wave  normal-shock  combustion 
tube  method,  it  does  not  appear  to  be  a 
method  for  determining  chemical  reaction 
kinetic  rates  for  individual  reactions.  The 
interactions  between  reactions  and  rates 
for  a  multicomponent  reaction  are  too  com¬ 
plicated.  However,  the  method  may  have  a 
distinct  advantage  in  verifying  heat  release 
rates  obtained  from  combinations  of  re¬ 
actions  and  rates  that  were  individually 
observed  from  other  special  experiments,  or 
for  determining  heat  release  rates  for  com¬ 
plicated  fuels,  such  as  hydrocarbons,  or 
mixtures  of  hydrocarbons,  such  as  JP-4  ,  As 
mentioned,  a  preliminary  experiment  with 
methane  and  vitiated  air  combustion  demon¬ 
strated  that  the  same  type  of  pressure  pro¬ 
file  data  can  be  obtained  as  with  H2-air. 
Also,  the  method  can  be  used  with  milch 
richer  mixtures  than  the  very  dilute  con¬ 
centrations  examined  in  the  usual  shock 
tube  studies. 

Conclusions 

The  experimental  method  of  measuring 
pressure  profile  history  of  a  chemically 
reacting  high  speed  gas  stream  yields  data 
which  substantiate  recent  results  reported 
from  shock  tube  experiments  on  chemical  re¬ 
action  schemes  and  their  rates  for  the  H2- 
air  reaction.  In  addition,  a  mathematical¬ 
ly  simple  method,  based  on  the  partial 
equilibrium  concept,  was  devised  for  calcu¬ 
lating  the  chemical  reaction  history  for 
H2-air,  and  the  results  were  found  to  be  in 
close  approximation  to  the  more  complicated 
reaction  schemes.  The  implication  is  that 
this  partial  equilibrium  (PE)  approach  can 
be  used  in  lieu  of  the  full  equations  for 
the  H2-air  system. 

Experimental  phenomena  similar  to  those 
observed  for  the  H2~air  reaction  were  also 
observed  for  the  methane-air  combustion, 
where  the  air  supply  was  45-percent  viti¬ 
ated  during  heating  to  approximately  3100°R. 
Ignition  delay  was  calculated  to  be  4  to  5 
times  faster  than  measurements  from  other 
investigators.  Rather  large  concentrations 
of  free  radicals  from  the  combustion  heated 
air  supply  may  be  an  explanation  for  the 
differences.  Thus,  a  nonvitiated  air  sup¬ 
ply  is  highly  desirable  for  combustion 
processes  where  the  effects  of  free  radicals 
are  difficult  to  determine. 


The  close  agreement  oi  the  experimental 
data  and  calculated  data  ior  the  H2~air 
system  also  implies  that  the  experimental 
technique  will  prove  to  be  satisfactory 
for  the  determination  of  heat  release 
rates  in  other  fuel-oxidizer  systems.  At 
present,  it  appears  to  be  a  useful  steady- 
state  technique  which  is  capable  of  sup¬ 
plying  empirical  equations  for  the  overall 
reaction , 

Results  of  this  type  of  experimental 
and  analytical  approach  to  analyze  hydro¬ 
carbon  fuel  heat  release  should  be  appli¬ 
cable  to  the  design  of  combustors  for 
supersonic  combustion  ramjets.  Another 
possible  application  is  in  simplifying 
computations  involving  simultaneous  mixing 
and  combustion  or  in  simultaneous  combus¬ 
tion  in  a  method-of-characteristics  net. 

The  method  is  not  perfected,  in  that 
(1)  the  entering  gas  mixture  was  not  uni¬ 
form,  and  (2)  shock  oscillations  were 
encountered  at  some  conditions .  Refine¬ 
ments  could  be  made  with  further  efforts. 
However,  the  technique  appears  promising 
to  obtain  heat  release  rate  data  for  gas¬ 
eous  reactions  where  little  or  nothing  is 
known  of  the  reaction  kinetic  scheme  or 
the  reaction  rates. 
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Appendix 


wher  e 


Analysis  of  the  Effect  of 
Vail  Boundary  Layer  and  Heat  Transfer1 

The  general  "influence  coefficient" 

Eq .  (14)  for  uniform  flow  in  a  constant 
area  duct  is  assumed  to  be  applicable  to 
fully  developed  flow  where  heat  transfer 
and  friction  exist.  However,  in  the  ex¬ 
periment  usirg  a  normal  shock  at  the  en¬ 
trance  to  a  tube  and  unreacting  gas,  the 
boundary  layer  growth  is  very  similar  to 
that  on  a  flat  plate;  therefore,  the  flow 
m  the  tube  was  treated  as  if  boundary 
layer  displacement  thickness  formed  an 
effective  area  constriction  in  the  constant 
area  tube.  Under  these  conditions,  heat 
transfer  to  the  cool  walls  serves  to  reduce 
the  boundary  layer  thickness. 

Thus  the  interior  flow,  excluding  the 
displacement  boundary  layer,  can  be  treated 
as  isentropic  when  no  combustion  exists  . 
From  consideration  of  conservation  of  mass, 
the  effective  area  (Ae)  can  be  written  in 
terms  of  the  displacement  thickness,  6*: 

A 

=  l  -  25*/r  (1-1) 

A  - 
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p  U  LI  dx 
1  e  e'e 


and 


ji  l  7 


dy 


where  p  is  the  density,  u  the  axial  velocity, 
p.  the  viscosity,  x  the  distance  from  the 
leading  edge  of  the  tube,  y  the  distance 
measured  normal  to  the  tube  wall,  6*  the 
boundary  layer  thickness,  f (q)  the  Illing- 
worth-Levy  stream  function  variable,  and  the 
subscript  (e)  denotes  the  interface  between 
the  boundary  layer  and  the  isentropic  core 
region.  By  using  the  perfect  gas  law  and 
the  Crocco  relation,  it  can  be  shown  that 
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f(n6) 
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(f')2  dTj 
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where  A^  is  the  tube  cross  section  area  and 
r  is  the  radius . 

In  the  entrance  region,  the  axial  pres¬ 
sure  gradient  is  relatively  small  and  has 
negligible  effect  on  the  growth  of  the 
boundary  layer;  therefore,  the  boundary 
layer  can  be  analyzed  in  the  same  manner 
as  the  boundary  layer  on  a  flat  plate. 

If  the  laminar  boundary  layer  equations 
for  zero  pressure  gradient  and  Prandtl  num¬ 
ber  of  unity  are  transformed  using  the 
Illingworth-Levy  transformation2  varia¬ 
bles  and  the  viscosity  is  assumed  to  vary 
linearly  with  temperature,  the  momentum 
equation  will  be  recognized  as  the  incom¬ 
pressible  flat  plate  equation  of  Blasius 
(see  for  example  Dorrance2) .  The  solution 
to  the  energy  equation  will  be  Crocco's 
Integral.2  The  displacement  thickness,  by 
definition,  is 
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and  under  the  transformation,  it  becomes 
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where  Tw  is  the  wall  temperature  of  the  tube, 
y  the  ratio  of  specific  heats,  Me  the  Mach 
number,  and  f'  the  Illingworth-Levy  non- 
dimensional  axial  velocity, 


From  the  Blasius  solution, 
shown  that 

\  t  f(n5)  “  1-22 
and 

rn6 


it  can  be 
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For  low  subsonic  velocities,  Eq .  (1-5) 
becomes 


<T  = 


1.73 
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"^Equations  derived  by  C.  Willbanks. 

2 

Ref.  Schlicting,  H.,  "Boundary  Layer 
Theory,"  4th  Ed.,  McGraw-Hill,  1960, 
p.  169. 

3 

Ref.  Dorrance,  W.  H.,  "Viscous  Hyper¬ 
sonic  Flow,"  McGraw-Hill,  1962,  p.  30-33. 


and  assuming  pe,  u^,  and  Te  constant,  Eq. 
(1-6)  can  be  simplified  further  to 


1 , 73x  _w 
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where 


Re  = 
x 


Pluixi 


The  subscript  (i)  denotes  conditions  at  the 
entrance  to  the  tube.  The  assumption  of 
constant  pe>  ue,  and  Te  is  justifiable  here 
because  the  correction  term  for  boundary- 
layer  is  small  and  also  because  these  quan¬ 
tities  do  not  change  greatly  from  one  end 
of  the  tube  to  the  other. 

Combining  Eqs .  (1-1)  and  (1-7)  yields 


A. 

l 


1 


3.46 

r 


Re  1/2  T, 
x 


(1-8) 


This  area  ratio  can  be  used  in  conjunction 
with  the  isentropic  flow  relations  to  de¬ 
termine  the  pressure  drop  through  the  tube 
caused  by  friction  and  heat  transfer . 

Clearly,  a  similar  analysis  could  be 
made  for  the  turbulent  boundary  layer. 
Moreover,  the  case  of  a  chemically  reacting 
gas  could  be  treated  in  an  analogous  manner, 
although  the  analysis  would  be  complicated 
by  the  necessity  of  considering  the  heat 
release  distribution  in  the  tube. 
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The  design  of  a  high  speed  airbreathing  propulsion  system  (SCRAMJET) 
may  depend  on  more  precise  information  on  chemical  kinetics  of  the  com¬ 
bustion  processes  than  is  now  available,  except  for  the  very  simplest 
reactions.  An  experimental  method  was  devised  whereby  the  macroscopic 
chemical  reaction  history  can  be  traced  from  wall  pressure  measurements 
in  a  high  speed  reacting  gas  in  a  tube.  The  chemical  reaction  is 
started  by  means  of  a  normal-shock  compression  at  the  tube  entrance 
through  which  premixed  fuel  and  air  is  heated  to  combustion  temperature. 
The  experiments  were  conducted  with  H2  and  vitiated  air,  and  the  data 
were  compared  with  calculated  pressure  profiles  using  recent  reaction 
kinetic  information  obtained  from  shock  tubes .  Good  agreement  between 
experiment  and  calculation  was  found  when  the  flow  was  steady  and  no 
oscillations  were  present,  indicating  that  the  kinetics  which  govern  the 
chemical  reactions  in  shock  tube  experiments  can  be  reproduced  with  the 
standing-wave  experiment.  Wall  effects  are  evaluated.  In  addition,  a 
semi-empirical  analysis  was  made  of  the  ignition  delay  and  recombination, 
with  the  assumption  that  the  two-body  reactions  are  essentially  in  equi¬ 
librium  during  the  slower  recombination  period.  The  mathematics  of  this 
approach  is  simple  and  is  potentially  useful  with  machine  computations 
involving  chemistry  and  mixing,  or  chemistry  and  aerodynamics.  Finally, 
a  short  experiment  with  methane  combustion  in  the  tube  produced  similar 
phenomena,  indicating  that  the  method  may  be  useful  for  determining 
rates  of  heat  release  for  gases  where  little  is  known  of  the  kinetics. 
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